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We reports direct and scalable measurement of multiparticle entanglement concun'ence and three-tangle with 
embedding photonic quantum simulators. In this embedding framework [Phys. Rev. Lett. Ill, 240502 (2013)], 
W-qubit entanglement monotone, which associates with non-Hermitian operators, can be efficiently measured 
with only 2 (for even N) and 6 (for odd N) local measurement settings. Our experiment uses a multiphoton 
quantum simulator to mimic the dynamical entanglement evolution and track its concurrence and three-tangle. 


Entanglement lies at the heart of quantum mechanics [1,2] 
and plays an important role in quantum information sci¬ 
ences including quantum simulation [3], quantum computa¬ 
tion [4, 5], quantum network [ 6 , 7], and quantum metrology 
[ 8 ]. Immense experimental advance has allowed the creation 
of growing number of entangled particles in laboratories. To 
date, entanglement among fourteen trapped ions [9] , eight 
photons [ 10 ] , and five superconducting qubits [ 11 ] have been 
demonstrated. 

Not only the generation of multiparticle entanglement, but 
also their verification and characterization, become exceed¬ 
ingly difficult for large number of quantum bits (qubits). 
Tracking the presence of multipartite entanglement has been 
considered as an important criteria for confirming quantum¬ 
ness in quantum machines [12]. A fundamental obstacle for 
quantifying entanglement [13] is that its function usually in¬ 
volves a complex-conjugation operator, which doesn’t asso¬ 
ciate with a Hermitian observable directly. 

The conventional method of quantum state tomography 
[14] requires ^4^ measurements for A^-qubit states, which 
is not scalable. Entanglement witness [15-17] has been de¬ 
signed to detect the presence of certain classes of multiparticle 
entanglement with polynomial complexity of measurement 
settings. However, it is not universal—as different classes of 
entangled states would need different witnesses—and doesn’t 
provide quantitative measures of entanglement. A single¬ 
shot measurement of concurrence [18] has been reported for 
a two-qubit state. Such methods, however, rely on collective 
measurements on multi-fold copy of the quantum states [19] 
which is experimentally demanding for multiparticle entan¬ 
glement. 

In this Letter, we report efficient measurement of con¬ 
currence [ 20 ] and three-tangle [ 21 ] entanglement monotones 
[ 22 ] evolution from initial separable states for two- and 
three-qubit system with embedded photonic quantum simula¬ 
tors. The photonic simulators are initialized in W- [23] and 
Greenberger-Horne-Zeilinger (GHZ) - [24] type entangled 
states, followed by local non-unitary operations, to mimic the 
dynamics of multiparticle entanglement. With only 2 (for 


even number of qubits) and 6 (for odd number of qubits) 
measurements, we reconstruct the entanglement monotone 
dynamics, confirming the high efficiency of the embedding 
framework. 

By definition, entanglement monotone [22] is zero for sep¬ 
arable states and doesn’t increase in average under local oper¬ 
ations and classical communication. The entanglement mono¬ 
tone for bi-particle system—concurrence—can be written as 
\{AoK) \ , where Aq = YY (here, X, Y, and Z are Pauli op¬ 
erators) and K is the complex-conjugation operator (K |'0) = 
|'0*)). The monotone for tri-particle system—three-tangle— 

is {AiKf + {A 2 Kf - {AsKf , where Ai = XYY, 
A 2 = ZYY, and A 3 = lYY (/ is the identify operator). 
These definitions can be generalized to multiparticle case [22]. 
Eor systems with even-iV particles, the monotone can be writ¬ 
ten in the form of concurrence, where Aq = Eor sys¬ 

tems of odd-A^ particles, the monotone is similar to three- 
tangle by replacing Ai, A 2 , and A 3 with Ai = XY^^~^, 
A 2 = ZY®^~^, and A 3 = It can be noted that 

the entanglement monotone of multiparticle system can be ex¬ 
pressed in the term of (AK). Generally, (AK) is not a real 
number as AK is not a Hermitian observable. 

Candia et al. [25] proposed a protocol for efficient mea¬ 
surement of multipartite entanglement with embedding quan¬ 
tum simulators (see Eig. 1). In this embedding framework, the 
wavefunction of multiparticle system IQ) is separated into real 
and imagine components and entangled with an additional an¬ 


cillary qubit, resulting in an enlarged state 


tp) as 


= |0) 0 Re IV') + |1) ® Im |i/;). (1) 

The real and imaginary parts of the value of (AK) can be 
directly measured by a Hermitian observable on the enlarged 
state by noting that 

Re {AK) = (ZA ), (2) 

Im (AK) = - (XA ), (3) 

where Z and X are the Pauli operators on the ancillary qubit. 
Therefore, for an even-Ai particle system it will take only two 






2 


(4'(f)M|4'*(0) = ? 

'F(0)) = Re|'F(0)) + Im|‘l'(0)) N qubits—^ - H = A + iB - ' ' l / 

' ' ' full tomography 


embedding 

embedding 

efficient 

coding 

simulation 

measurement 


ancillary qubit — 

TO,)H0)R.|T,0)) + |1)I.|T,0|) 


- f iB 

iA\ 

H = \ 



iB) ■ 


('I'(0|Z^|'RW)-/('I'(0|X4|'I'(0) 


FIG. 1; Embedding simulator for efficient entanglement measure. The entanglement in original dynamic system H cannot be directly measured 
as it contains complex terms like (fi/| A |T'*). Through embedding the initial state and Hamiltonian into an enlarged simulator assisted by an 
ancillary qubit, the complex terms are allowed to efficiently measure on the enlarged state. 


Local operations Measurements 



FIG. 2: Experimental setup. Ultraviolet (UV) femtosecond laser 
pulses (394 nm, 120 fs, 76 MHz) pass through two type-II BBO crys¬ 
tals to produce two pairs of entangled photons. Two single photons, 
one from each pairs, are mixed on a PDBS to generate W-type entan¬ 
gled state. To simulate the three-tangle system, the PDBS is switched 
to a PBS to generate GHZ-type entangled four photons. The multi¬ 
photon sources are sent to next stage to perform local operations by 
HWPs and a polarization dependent Mach-Zehnder interferometer. 
Finally, the photons are measured in by four-fold coincidence count¬ 
ing. BBO, beta-barium borate; PDBS, polarization dependent beam 
splitter with a transmission of 0.72 for H photons and 0.28 for V 
photons ; PBS, polarization beam splitter; HWP, half-wave plate; 
QWP, quarter-wave plate. H and V denote horizontal and vertical 
polarization. 


measurements of observables ZAq and XAq on the enlarged 
state to evaluate the entanglement, and for an odd-A^ particle 
system it needs only six measurements of observables ZAi, 
XAi, ZA 2 , XA 2 , ZAa, and XA 3 . 

By embedding the initial wavefunction |'0(O)) and quantum 


dynamics H into an enlarged simulator initialized in 


m 


with an appropriate Hamiltonian H, the complex-conjugation 
expectation relation at various time still holds as 


{AK){t) = {ZA){t)-i{XA){i), (4) 


where 

" = + is) ® 


with real matrix A'^ = A and = —B. This allow us to 
track the multiparticle entanglement dynamics efficiently. 


Now we proceed with experimental demonstration using 
multi-photon quantum simulators [26]. To demonstrate the 
general working principle for multiparticle systems of even 
and odd number of qubits, we choose to implement the mea¬ 
surement of dynamical concurrence of a two-qubit entangled 
state and three-tangle of a three-qubit entangled state. The 
first example, a concurrence system, starts from an initial state 
|'0(O)) = |0) 0 |0) with a Hamiltonian H = 2f0F + X 0 0. 
The task is to evaluate the concurrence at arbitrary time t as¬ 
sociated with the unknown quantum state |V'(f))- According 
to embedding framework, we set an enlarged embedding sim¬ 


ulator with an initial state 


' 0 (O)y = | 0 ) 0 | 0 ) 0 | 0 ) and a 


Hamiltonian H = I ® X ®Y — Y ® X ® Z. 


We implement this enlarged dynamical system in a three- 
photon compiled simulator. Our experimental setup is showed 
in Fig. 2. We use single photons produced from spontaneous 
parametric down-conversion [27] as qubits, where the hori¬ 
zontal {H) and vertical (V) polarization are used to encode 
|0) and 11), respectively. The compiled simulator runs in three 
stages; ( 1 ) create an initial entangled source, ( 2 ) perform local 
operations, (3) finally, readout the expectation of observables 
XAq and ZAq. 

At stage (1), three single photons (ancillary, 1, and 2) are 
entangled in W-type state (|011) — |101) + |110))/-\/^ by 
post-selecting the photon 3 in |0) after mixing two pairs of 
Bell-state (|01) + |10))/-\/2 photons on a polarization de¬ 
pendent beam splitter (PDBS) with a transmission of 0.72 
for H photons and 0.28 for V photons [28]. At stage (2), 
the entangled photons are sent to perform local operations by 
half wave plates (HWPs) and a polarization dependent Mach- 
Zehner (MZ) interferometer. The ancillary photon was oper¬ 
ated with X operation and the photon 2 was operated with 
Z and X operations by passing through HWPs. Now the 
photons state are (|000) + |011) — |110))/-\/3. Then pho¬ 
ton 1 passes through a polarization dependent MZ interfer- 

, where the 

coming photon is split into two spatial modes by a polariza- 


ometer operated with 


cos y/^t 0 
0 Smy/2t/y/2 
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FIG. 3: Experimental results of entanglement dynamics in concurrence system. (a)(b) The time-dependent expectation of observable ZYY 
and XYY. (c) The concurrence computing from two measurements of ZYY and XYY. The error bars represent poissonian statistics in 
photon counting. The red curves are theoretical simulation counting for the imperfect interferometric visibility. 
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FIG. 4: Experimental results of entanglement dynamics in three-tangle system, (a)-(f) The time-dependent expectation of observable ZXYY, 
XXYY, ZZYY, XZYY, ZIYY, and XIYY. (g) The three-tangle computing from previous six measurements. The error bars represent 
poissonian statistics in photon counting. The red curves theoretical simulation counting for the imperfect interferometric visibility. 


tion beam splitter (PBS) regard to its polarization, and each 
pass through a HWP, then combine on a second PBS be¬ 
fore anther HWP at the output port. After this stage, it will 


generate the desired enlarged state 


'0(f)), in the form of 


cos(V2f) |000) -F sin(V2f)/V2 |011) - sin(y2f)/y2 |110). 


We carry out the measurements on stage (3). By succes¬ 
sive stepping the orientation of HWPs inside the MZ inter¬ 
ferometer, the time parameter f-dependent evolution is set at 
the points; 0, Af, A2f, A3f, •••, where Af = 7r/12-\/2- 
To evaluate the concurrence and reconstruct its dynamics, 
we made two measurements of ZYY and XYY at various 
time t. The experimental results are shown in Fig. 3(a) and 
(b). Then we computed the concurrence directly and shown 
it in Fig. 3(c). The concurrence dynamics is well fitted by 
\a ■ sin(-\/2f)/'\/2| with a = 0.59 ± 0.01, agreeing with the 
theoretical simulation scaled by an amplitude factor a count¬ 
ing for imperfect photonic interference visibilities in the cre¬ 
ation of three-photon W state [28]. 


Next we present the second example, a three-tangle system 
with an initial state of |0(O)) = |0) ® |0) ® |0) and a Hamil¬ 
tonian of H = X ® X Z) X. To measure the entanglement 


of |'0(f)), the enlarged state '0(f)y are implemented with an 
embedding simulator of initial state 0(0)^ = |0000) with 
Hamiltonian H = —YXXX . 


The experiment was also run in three stages. At stage (1), 
the initial entangled source (|0110) + |1001))/'\/2 is created 
by combining two pairs of Bell-state (|01)-F|10))/ photons 
on a PBS. Then at stage (2), the photon 1 and 2 pass through 
a HWP to implement an X operation. The photon 3 passes 
through a HWP to implement a Z operation. And further, a 

f cosf 0 \ , 

non-umtary operator I ^ s' t I applied to photon 1 

by passing through a polarization dependent MZ interferom¬ 
eter. These two stages will make the desired enlarged state 
'0(f)^ = cos(f) |0000) — sin(f) At the final stage 

(3), we measured six observables ZXYY, XXYY, ZZYY, 
XZYY, ZIYY, and XIYY at time 0, At, A2f, A3f, • • •, 
where At = 7r/12. The Fig4 (a)-(f) shows the experimen¬ 
tal results. So we can directly compute the three-tangle. Fig. 
4(g) show the evolution of three-tangle, theoretically fitted by 
\a ■ sin(2f)|^ with a = 0.70 ± 0.01, which is the visibility 
of the generated four-photon GHZ state in the superposition 






















4 


basis. 

These two examples showed that the embedding framework 
can be used for efficient measurement of multiparticle entan¬ 
glement dynamics. We have observed the dynamical entan¬ 
glement oscillation well according with the theoretical pre¬ 
diction. The experimental imperfect performance of reduc¬ 
ing the concurrence oscillation amplitude mainly comes from 
the high order emission of spontaneous parametric down- 
conversion photon pairs [2]. 

In conclusion, we have presented a proof-of-principle 
demonstration of efficiently measuring entanglement dynam¬ 
ics for concurrence and three-tangle dynamics system on com¬ 
plied embedding photonic simulators. We have proved the 
high efficiency and feasibility of embedding framework for 
multipartite entanglement measurement. This paradigm can 
be applied to other physical systems such as trapped ions [9] 
and superconducting qubits [11] for scalable entanglement 
measurement, and may offer a promising way to study the en¬ 
tanglement of large-scale dynamical multipartite system. The 
embedding framework, exploiting additional qubits to reform 
the original wavefunction information, may inspire other ex¬ 
citing applications, such as simulation of non-unitary Majo- 
rana dynamics [29]. 

Note: When preparing our manuscript, we became aware 
of a related experiment work from Andrew White’s group. 
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